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We report a laboratory investigation of weak turbulence of water surface waves in the gravity¬ 
capillary crossover. By using time-space resolved profilometry and a bicoherence analysis, we observe 
that the nonlinear processes involve 3-wave resonant interactions. By studying the solutions of the 
resonance conditions we show that the nonlinear interaction is dominantly ID and involves collinear 
wave vectors. Furthermore taking into account the spectral widening due to weak nonlinearity 
explains that nonlocal interactions are possible between a gravity wave and high frequency capillary 
ones. We observe also that nonlinear 3-wave coupling is possible among gravity waves and we raise 
the question of the relevance of this mechanism for oceanic waves. 


A large ensemble of nonlinear waves can exchange en¬ 
ergy and develop a turbulent state. The statistic prop¬ 
erties of such wave turbulence have been described the¬ 
oretically for weak nonlinearity in the framework of the 
Weak Turbulence Theory (WTT). In this theory, only 
resonant waves are able to exchange significant amounts 
of energy over long times due to the weak nonlinear cou¬ 
pling. The predicted phenomenology of the stationary 
statistical states resembles that of fluid turbulence: en¬ 
ergy is injected at large scales and cascades down scale 
to wavelengths at which dissipation takes over and ab¬ 
sorbs energy into heat. A major difference with fluid 
turbulence is that analytical predictions for the station¬ 
ary spectra (and other statistical ^antities) can be de¬ 
rived for weak wave turbulence [ll-l3|. Sea surface waves 
are among the pioneering plwsical systems that led to the 
development of the theory . The theory was applied to 
a vast amount of waves (in plasmas [l|, solar winds [H, 
nonlinear optics ^antum superfluid vortices 0 , vi¬ 
brated elastic plates | 8 |,...). 

For isotropic systems, the predicted energy spectrum 
E(k) has the following expression 


E{k) = 


( 1 ) 


where /c = |k| is the wavenumber, P the energy flux, 
C a dimensional constant that can be calculated and a 
the spectral exponent. N is the number of waves tak¬ 
ing part in the resonances. Usually N — 1 corresponds 
to the order of the nonlinear coupling term of the wave 
equation [N = 3 for quadratic nonlinearities, = 4 for 
cubic ones,...). The waves have then to satisfy the reso¬ 
nance conditions such as ki = k 2 + ks and cji = co ’2 + (^3 
(for 3-wave interaction). In some cases, these resonances 
conditions do not have solutions. This is the case in par¬ 
ticular for gravity waves at the surface of water. The 
dispersion relation (for infinite depth) is cj = and 

its negative curvature does not allow for solutions of the 
resonance conditions for 3 waves. Thus the resonances 
are expected to involve 4 waves. At small wavelengths, 
water waves are capillary waves for which the dispersion 
/ \ 1/2 

relation is cj = f whose curvature allows for 


3-wave resonances (7 is the surface tension and p the 
density). The predicted spectra for water waves are thus 
expected to be E{k) oc for capillary waves 

and EiP) oc Q for gravity waves. Labora¬ 

tory experiments largely fail to reproduce this predic¬ 
tions, in particular for the gravity waves. In large or 
small waves tanks, the spectral exponent of the gravity 
waves is seen to vary strongly with the forcing intensity 
and to be close to the WTT predictions at the highest 
forcing mamitude, at odds with the weak nonlinearity 
hypothesis Furthermore the measurement of the 

exponent of the injected power is also different from that 
predicted by WTT 0 - Recent work on water waves 
and vibrated plates suggest that wideband dissipation is 
most likely responsible for the latter observation [uHIi- 
Another experiment also suggest that several regimes of 
wave turbulence of water wave may exist depending on 
the intensity and frequency of the forcing 0 - Neverthe¬ 
less the question of the order of the interaction remains 
of prime importance to test the theory. In this letter, 
we report a high order statistical analysis that directly 
probe the nonlinear interaction among waves. We im¬ 
plement a time resolved, 2D profilometry of the water 
surface deformation. The accessible wave lengths cor¬ 
responds to capillary waves and to the gravity-capillary 
crossover at which the order of interaction supposedly 
switches from 3-wave at small wavelength to 4-wave at 
large wave lengths. We investigate to resonant interac¬ 
tion and the impact of the wave amplitude on the energy 
transfers. 

The experimental setup consist of a rectangular plas¬ 
tic vessel of 57 x 37 cm^ filled with 10 1 of water to a 
rest height ho = 5 cm (fig. [T]). Two curved walls are 
used to improve the isotropy by divergent waves reflex¬ 
ion. Surfaces waves are excited by horizontally vibrating 
the vessel at frequencies in the range [0.5,1.5] Hz. Waves 
are measured using the Fourier transform profilometry 
technique which enables a full space-time characteriza¬ 
tion of the waves [l5|. We used filtered softened water 
and thoroughly washed the tank. Similarly to Przadka et 
al. 0 we use anatase titanium dioxide particles (Kro- 
nos 1001 ) that do not alter the measured water surface 
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FIG. 1. Sketch of the experiment (Top view). The vessel is 
translating horizontally with a constant amplitude frequency 
modulated oscillation, which frequency is chosen in the inter¬ 
val [0.5,1.5] Hz. Two curved walls are placed on left and right 
sides to improve the isotropy by diverging waves reflexion 


tension and do not induce additional dissipation at the 
surface. Thanks to this pigment, it is possible to project 
a pattern at the very surface of water. When the wa¬ 
ter surface is deformed, the pattern seen by a camera 
is changed. The alteration of the pattern can then be 
inverted to recover the deformation of the surface El- 
Here, the deformation of the pattern is recorded by a 
high-speed camera over a 20 x 20 cm^ surface at the cen¬ 
ter of the tank, with 1024 x 1024 pixels resolution at 
250 frames/s. Datasets are made of 15 movies which du¬ 
ration is 87 s for each movie. 

We show in fig. [2] the space-time power spectrum 
of the velocity field v = ^ where rj{x,y,t) is 
the altitude of the water surface. We first compute 

Eyk,c.) = (Kk,c.)|2) (2) 

where i;(k, cj) is the space and time Fourier transform of 
the velocity. The time Fourier transform is computed by 
selecting a 16 s time window of the signal. The average 
(...) is a time average over the time windows. oo) is 

then summed over directions of the 2D wavevector k to 
provide a 2D picture of E'^{k,uj) (fig. [2]). Energy is seen 
to be concentrated around the linear dispersion relation 
of gravity-capillary waves 

/ \ 1/2 

uu = (gk -h —k^j (3) 

The isotropy of E^(k, cj) is shown in the inset of fig. [2] 
at a given frequency of 10 Hz. Energy is convincingly 
distributed over all directions. The energy concentra¬ 
tion around the dispersion relation is due to nonlinear 
spectral widening predicted by the WTT framework. 
Note that no secondary branches of the dispersion re¬ 
lation are seen contrary to what was reported in Hi- 
Our regime corresponds to the second regime of turbu¬ 
lence reported in [1^ at the weakest magnitude of the 
waves. The wave steepness of our data is indeed small: 

a = ^ Js II V/i(x, t)||^ dxdy^ = 0.025 thus our wave 

field is weakly nonlinear. 
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FIG. 2. Space-time Fourier spectrum of the velocity field of 
the waves E'^{k^uj) (see text for definition). The color scale 
is log^Q E^(/c, cj). The solid black line is the theoretical deep 
water linear dispersion relation for pure water = gk-\- ^k^ 
with 7 = 72 mN/m. Energy is localized on the dispersion 
relation and can be observed for frequencies up to 60 Hz. 
The cr ossove r between gravity and capillary waves occurs at 
= \/~pg~hi — 1207r corresponding to a wavelength of 1.7 cm 
and a frequency equal to 13 Hz. Inset: E(k,ct;) at Ct;/27r = 
10 Hz. The energy distribution is fairly isotropic. The black 
circle corresponds to the linear dispersion relation. Golor: 
logj^Q E(k, ct;) between —8 and —4. 


As discussed in the introduction, the order of nonlin¬ 
ear wave interaction depends on the order of the nonlin¬ 
earity and the possibility or not to have solutions for 
the resonance equations. McGoldrick and then Sim¬ 
mons EiiDi has investigated the 3-waves resonant solu¬ 
tions in gravity-capillary regime. A 3D representation of 
the solutions for a given wavevector k 2 is shown in fig. [3l 
As the curvature of the dispersion relation changes sign, 
solutions of the 3-wave resonance conditions exist even 
if the wave vector k 2 lies in the gravity range. Such a 
solution exist only if the norm of the vector ks exceeds 
a minimum value as can be seen from the red line in 
fig. [3l This minimum value of /cs is reached when the 
three wave vectors ki, k 2 and ks are collinear (ID sit¬ 
uation). This can be translated in the frequency space: 
at a given frequency Cc; 2 , there is a minimum value of 
Co’s that allows for a resonant wave at uji. In the fol¬ 
lowing, we note this minimum value (fig. [31(b)) 

and co’f^*’^(co’ 2 ) the corresponding value of uji (fig. |3l(c)). 
These values of the resonant frequencies correspond to ID 
wave interaction. By scanning values of co’ 2 , one observes 
that there is a overall minimum value for resonant 
corresponding to the degenerate case of Wilton wave : 
uji = 2uj2 = 2cc;s = 27r x 19.6 Hz (fig.[3lfc) ) 1^ . Measure¬ 
ments performed by Hammack et al. [22| confirm the 
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FIG. 4. Representation of 3-wave coherence \C{uJi^U 2 ^o^z)\ 
(see text for definitions). The specular light blue pattern on 
the background corresponds to statistical convergence noise so 
that the statistical convergence is about 10“^’^. Here uj/2 ti is 
chosen equal to 10 Hz. Thus the observed line of correlation 
lies above the convergence level by more than one order of 
magnitude at the maximum of correlation. This is the sign 
of the presence of significant 3-wave coupling in the signal. 
The line lies on the resonance curve uji = 002 + (^3 as expected 
from such statistical estimator. Color is log^^Q \C\. 


FIG. 3. Solutions for the 3 wave resonances of gravity capil¬ 
lary waves, (a) The red surface corresponds to the dispersion 
relation (|3]). First a given value of k 2 is chosen. The blue sur¬ 
face corresponds then to the sum Ct;(k 2 ) +Ct;(k 3 ). The resonant 
wave at ki must verify the resonance conditions ki = k 2 + ka 
and (jUi = (jU 2 (^3 and thus lie on the intersection of the red 
and blue surfaces (thick black line). The red line at the bot¬ 
tom of the axes corresponds to the projection of the black 
line, (b) Evolution of uj^^'^(uj 2 ) that is the minimum value of 
(jU 3 for which there is a possible solution for a given value of 
UJ 2 . (c) corresponding value . 


existence of these 3-waves resonant coupling. 

To investigate the 3-wave coupling in our experiment, 
we study third order correlations of the velocity field. 
From v{x, y, t), we compute the Fourier transform in time 
over 4 s time windows so that to obtain v{x^y^uj). Cor¬ 
relations are then computed as 


C(ci;i, CJ 2 , ^3) 


y, uii)v{x, y, ui 2 )v{x, y, 0 ) 3 )))| 
[E-{LUi)E-{uj2)E-{LU3)f^ 


( 4 ) 

where * stands for complex conjugation and the aver¬ 
age ((...)) stands for an average over the time windows 
and a space average over (x, y) positions on the im¬ 
age. = ((|'i;(x, cj)!^)) is the frequency spectrum. 


With such a normalization, the coherence lies between 0 
(no correlation) and 1 (perfect correlation). For a sta¬ 
tionary (in time) signal such a third order correlation 
is expected to be nonzero only along the resonance line 
uji = co ’2 + <^ 3 - In order to check whether 3-wave cor¬ 
relation is truly observed, we show the coherence map 
at a given frequency (co’2/27r = 10 Hz) in fig. [H A line 
of correlation emerges from the statistical convergence 
noise that confirms that 3-wave resonant processes are 
indeed present in the signal. Such a line can be observed 
at all values of co’ 2 - Note that the correlation level is 
not homogeneous along the line showing preferential in¬ 
teraction among the waves. In order to look deeper in 
the nonlinear dynamics, we focus in the following on the 
bicoherence, defined as 


B{uJ2,UJ3) = C{UJ2 -\-UJ3,UJ2,UJ3) ( 5 ) 

It corresponds to the extraction of the coherence observed 
on the resonant line of fig. [4] which has been checked to 
be significant and above the convergence noise level. 

Figure [5] displays the bicoherence map corresping to 
the extraction of the coherence for all values of ci; 2 . A 
noticeable organisation of the coherence is clearly visible. 
First a wide curved line of high coherence (about 10“^) is 
observed. The crest of coherence lies on the curve 
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FIG. 5. Bicoherence B{uj 2 ^ujz) = C(ct ;2 + Ct; 2 ,c^a). The 

solid black line corresponds to the minimum value 
allowing exact resonances (see text). The dashed line corre¬ 
sponds to the transformation of the solid line when allowing 
for a (5A:/27r = 5 m“^ uncertainty on the resonance condition 
(see text). 

(solid black line) corresponding to unidirectional interac¬ 
tion of the waves as explained above. This observation 
suggests that the core of the nonlinear interaction in our 
experiment is quasi ID. The isotropy of the space-time 
spectrum of fig. 1 may be somewhat artificial and mostly 
due to chaotic mixing of the waves due to the boundary 
conditions (reflexion of curved walls) rather than nonlin¬ 
ear directional spreading. 

Another region of high coherence is seen also for weak 
values of and large values of UJ 2 (bottom-right part of 
fig. [5] along the axis). This region is made of two hori¬ 
zontal lines of high correlation values (a symmetric region 
exists next to the vertical axis). The discrete character 
of these two lines suggests that it corresponds to interac¬ 
tions of capillary waves with the first modes of our finite 
domain. Note that the eigenfrequencies of these modes 
lie in the gravity domain (0.9 Hz, & 1.7 Hz). The ex¬ 
istence of these regions seems paradoxical as no exact 
resonant interaction of waves is possible below the curve 

, ,min 
CJ3 . 

To explain this paradox, let us note that the coher¬ 
ence line along has a finite (nonzero) width. This 

width may be associated to the similar width of the en¬ 
ergy concentration along the dispersion relation of fig. O 
This width is due to nonlinear effects: indeed the coher¬ 
ence of the linear waves is altered by the nonlinear en¬ 
ergy exchanges with the other waves [i, [m . The width 
of the dispersion relation is thus a measure of the degree 
of nonlinearity of the system. Let us assume that 5k is 
the width in wavenumber of the dispersion relation {5k 
is taken constant for simplicity). It can be interpreted 


as a nonlinear uncertainty on the determination of the 
wavenumber. Thus the resonance condition translates 
into an inequality |ki — k 2 — ksj < 5k rather than an 
equality to zero. For simplicity, let us assume that as 
suggested by the data in fig. [5] the nonlinear coupling oc¬ 
curs solely through a ID mechanism. We look for the 
solutions = /c 2 + ± 5k and uji = 6^2 + 0 ^ 3 . The two 

corresponding solutions are shown as the dashed lines in 
fig.[5]for 5k/27r = b m“^. This value of 5k is reasonable in 
view of the width of the dispersion relation in fig. O The 
two dashed lines spectacularly encircle the blue regions 
of very low coherence and highlight the regions of high 
coherence. The new region of high coherence delimited 
by the dashed lines incorporates all previously discussed 
regions. This observation removes the apparent para¬ 
dox. When one takes into account the nonlinear spectral 
widening, the possible interactions are much more nu¬ 
merous, including frequencies well below the curve. 
In particular it opens a region of strongly nonlocal inter¬ 
action of a very low frequency gravity mode and 2 much 
higher capillary waves. For example, the mode at 0.9 Hz 
can interact with the whole interval of frequencies and 
thus initiate the energy cascade. 

Furthermore, as can be seen in the bottom left cor¬ 
ner of the picture, very low frequency modes can also 
interact among each other. As discussed in the intro¬ 
duction, the nonlinear interaction of gravity modes is 
usually assumed to involve 4 waves but incorporating a 
small nonlinear spectral widening may allow for 3-wave 
quasi resonances that may actually be more efficient to 
transfer energy. This possible interaction remains to be 
quantified precisely as in our simplified calculation, the 
widening is taken constant whereas it depends usually 
on the frequency If confirmed at larger scale, this 
mechanism could also be responsible for the discrepancy 
of the observed spectral exponents in the gravity range 
as compared to the WTT predictions that assume 4-wave 
interactions. 

In conclusion, we have shown that weak turbulence of 
water waves near the gravity-capillary crossover relies on 
a quasi ID 3-wave interaction. The theoretical reason 
for the selection of this ID mechanism while 2D reso¬ 
nances are a priori allowed remains to be investigated. 
Taking into account the nonlinear spectral widening into 
the study of the resonance conditions appear to be of 
prime importance at it changes significantly the range of 
possible interactions. In particular it allows for strongly 
nonlocal interactions between gravity and capillary waves 
and for 3-wave interaction among gravity waves, previ¬ 
ously assumed to be impossible. 
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